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Abstract 

We study the application of the joint resummation formahsm to Higgs production via gluon- 
gluon fusion at the LHC, defining inverse transforms by analytic continuation. We work 
at next-to-leading logarithmic accuracy We find that at low Qt the resummed Higgs Qt 
distributions are comparable in the joint and pure-Qr formalisms, with relatively small 
influence from threshold enhancement in this range. We find a modest (about ten percent) 
decrease in the inclusive cross section, relative to pure threshold resummation. 
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1. Introduction 



The gluon-gluon fusion process gg — > hX is expected to be the dominant mechanism for the 
production of rather hght Higgs bosons at the LHC [1]. Next-to-leading order (NLO) QCD 

corrections to the total production cross section for this process are known to be large, of the 
order of 70% [2]. This made the computation of NNLO corrections to the process an important 
task, which was completed very recently [3]. The NNLO corrections are found to be substantial 
as well, even at LHC energies, where they increase the total cross section by an additional 30% [3] 
over NLO. The dominant part of these contributions is directly related to soft and colhnear gluon 
emission, as demonstrated in Rcfs. [4, 5]. 

Soft gluon emission is manifested through the presence of logarithmic corrections in expres- 
sions for partonic subprocesses. These corrections arise from cancellations between real and 
virtual contributions to the differential cross section at each order in perturbation theory. They 
can be large when a cross section is sensitive to their behavior near the boundary of phase space. 
This is the case for threshold corrections where terms of the form a" ln^"~^(l — z)/{l — z) grow 
large in the limit z = Q'^/s 1, that is, when the partonic center-of-mass energy s approaches 
the invariant mass Q of the Higgs particle. Similarly, recoil corrections appearing in transverse 
momentum distributions, Q;"ln^"^^((5^/(5T)/(5r, become large when the Higgs boson is produced 
with a small transverse momentum Qt <^ Q- Because both the total production rate and the 
production characteristics, inchiding the transverse momentum distribution, will influence the 
search and analysis strategies for the Higgs, a close study of higher-order corrections in the soft 
and/or collinear limit is desirable. 

Resummation techniques for Higgs production have been established separately in the thresh- 
old [6] and in the recoil [7, 8, 9] cases, based on earlier work [10, 11, 12] on the related case of 
Drell-Yan dimuon production. A new analysis for threshold resummation, including the re- 
cently calculated [4] next-to-next-to-leading logarithmic (NNLL) coefficients C*^^^ and D'^'^\ can 
be found in [13, 14]. The resummcd Higgs-Qr distribution has also recently been reinvestigated 
in [15, 16, 17], using the NNLL B^'^^ coefficient calculated in [18, 19], and in [17] matched to the 
Higgs cross section at fixed Qt at NLO [20]. A complete understanding of soft gluon effects in 
differential distributions also requires a study of the relation between the two sets of corrections. 
A joint treatment of these corrections was proposed in [21, 22]. It relies on a refactorization 
of short-distance and long-distance physics at fixed transverse momentum and energy [22]. In 
this formalism, resummation of logarithmic corrections takes place in the impact parameter b 
space [12], Fourier conjugate to transverse momentum Qt space, and at the same time in Mellin- 
N space [10, 11], conjugate to z space. This guarantees simultaneous conservation of energy and 
transverse momentum of the soft radiation. At present, the joint resummation formalism has 
been developed to next-to-leading logarithmic (NLL) accuracy. A full phenomcnological study of 
the joint resummation formalism, as applied to electroweak Z production, was undertaken in [23]. 

It may seem surprising that threshold resummation is relevant in Higgs production at LHC 
energies, given that partonic threshold for the production of a Higgs at rest involves partonic 
fractions of the order of 10~^, even smaller than for the Z boson at the Tevatron, where the 
numerical effect of threshold resummation is rather small [23]. Nevertheless, the results of [3, 4, 5] 
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and [6, 13, 14] show that this is the case. Much of the difference between the Higgs and the Z 
boson can be attributed to the difference between the color charge Ca for the gluons that initiate 
the former and Cp for the quarks that initiate the latter. This difference reflects itself in an 
increased sensitivity to Sudakov logarithms associated with partonic threshold for gluon-induced 
processes. Gluon-induced processes, however, are also potentially sensitive to radiation at low 
X. Indeed, the steep x-dependence of the gluon distribution functions at low x enhances their 
sensitivity to threshold singularities. By examining the Higgs cross section at measured Qt, we 
will find an intriguing interplay between these two regions, which should have important influence 
of the shape of the resummed distribution. 

In any implementation of a resummation formula one has to address the issue that the run- 
ning strong coupling in the resummed expression is probed at virtualities reaching down to Aqcd 
and below. A number of methods have been proposed to deal with this singularity, while leaving 
intact the full resummed cross section [24, 25], relying on the analytic structure of the resummed 
cross section. In our previous study of Z boson production, [23], we made use of the minimal 
prescription [25] for the Mellin-inverse, and the prescription of [26] for the inverse Fourier trans- 
form, diverting it into the complex 6-plane according to the phase structure of the integrand (see 
below). In this way, we obtained a well-defined resummed cross section, valid for all nonzero 
Qt, which described the Qt distribution of Z bosons produced at the Tevatron. A complete fit 
requires a fitted non-perturbative parameter in the form of a Gaussian in b space [9] . In general, 
however, the effect of the nonperturbative parameter is quite modest, as also found in Ref. [27]. 
Toward large Qt, a matching procedure between the fixed-order and the resummed results led to 
good agreement with the data there as well. The treatment of the inverse 6-space transform can 
be carried out as well in pure Qt resummation as in joint resummation. We can therefore com- 
pare the two, to better understand the relationship between threshold and transverse momentum 
resummation. 

In this paper we study Higgs production via gluon-gluon fusion for both Qt and joint re- 
summation. The leading order cross section is of order and is mediated via a heavy quark 
loop. Because the Higgs couples to fermions proportionally to their masses, we make the by-now 
standard simplification of keeping only top quark loop contributions to the process. We also 
treat the mass of the Higgs boson as much lighter than twice the mass of the top quark. Under 
these conditions, the Higgs-gluon interaction can be described in the language of an effective 
Lagrangian [2] , replacing the top quark loop by a simple effective ggh vertex. It has been shown 
that the heavy top quark approximation is accurate within a few percent in the case of NLO cal- 
culations [6] for actually a much wider range of Higgs masses. Owing to similarities between the 
Drell-Yan production mechanism and the gluon fusion mechanism for Higgs production (under 
the above assumptions), the joint formalism developed in [23] can be fairly easily translated to 
the latter case. It turns out, however, that there are additional effects in the Higgs resummation 
case that require a more detailed analysis, as we will show. As suggested above, they are related 
to the gluon initial state for the reaction, resulting in amplified resummation effects through the 
basic replacement Cp — 4/3 — > = 3 in the step from Drell-Yan to Higgs, and to the important 
role of 1/x terms in the splitting functions. 

The joint resummation formalism for Higgs production in N and b space is described in 
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section 2, along with its relations to threshold and Qt resummations. In section 3, the method 
for inverting transforms is reviewed, along with the relation to more standard Qt resummations. 
Section 4 presents our numerical results for the joint and pure Qt cases. In section 5 we discuss 
the total (i.e., QT-integrated) cross section and the possible role of terms not brought under 
control by joint resummation. Finally, we draw conclusions. 

2. Joint resummation for Higgs production 
2.1. The joint ly-resummed cross section 

We start from the general formula for joint resummation [22, 23], applied to Higgs production: 



dQ^d?QT ^'^^''Jc^2m J {271 



X Cg/A{Q, 6, N, II, lip) exp [^J/(iV, 6, g, /x) ] Cg/siQ: b, TV, /x, fxp) ■ (1) 

This is the same general form as for vector boson production, but adapted to the Born cross 
section for Higgs production in the heavy top mass limit, rrih <C 2mt. The Higgs mass is rrih, and 
r = Q'^/S, where Q, which is set equal to m/j in the approximation associated with threshold and 
Qt resummation, sets the scale of the hard scattering. As is characteristic of joint resummation, 
the cross section is a double inverse transform, in impact parameter b and Mellin moment N. 
It resums logarithms in b and N to the level of NLL. Let us first summarize the content of the 
factorization implicit in this expression. 

Each of the factors in Eq. (1) is associated with a specific range of pcrturbative dynamics 
involved in the annihilation process. These ranges are set by the physical scales that enter the 
problem. Prom soft to hard these are, Aqcd at the scale of hadronic binding, followed by scales 
introduced by the transforms: Q/N and 1/6, followed by nih — Q, rrit and, finally s. We introduce 
a function x{bQ, N), to be specified in section 2.3., that interpolates between in the threshold 
limit, N ^ bQ and bQ in the limit of large impact parameter, bQ ^ A^, and is chosen in a 
way that makes it possible to generate logarithms in 1/6 and Q/N at the level of next-to-leading 
logarithm in each. At scales lower than Q/x, radiation is fully allowed, while beyond these scales, 
the transforms in N and b constrain final-state radiation. The resummations appropriate to these 
scales are specified in terms of finite-order coefficients at each logarithmic order. 

2.2. Coefficients and anomalous dimensions 

The shortest distances, with an off-shellness of order rrit, are represented by the perturbative 
function H{as). It collects the effects of hard virtual corrections and reads to first order [2, 18]: 

H{as) = 1 + ^ (27r' + n) • (2) 
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This is multiplied by the overall normalization constant 

^ 576^^ 

corresponding to the the heavy top mass limit, ruh <^ 2mt, with Gp the Fermi constant. 

At the longest scales, the functions C represent evolution between Aqcd and Q/x{bQ,N). 
They are given by 

Cg/H{Q, b, N, /X, fiF) = E Cg/j {N, asiQ/x)) {N, Q/x, /^f) fk/H{N, /x^) . (4) 

Here, the fj/H{N, /ip) are the parton distribution functions for hadron H at any convenient fac- 
torization scale /If- The Cg/j may be determined from the coefficients of the 5^{Qt) contributions 
in the order correction to the Higgs cross section [12], given our choice of if^^^ in Eq. (2). 
They read 

Cg,g{N,as) = l + ^7^^ (5) 
a 1 

Cg/,{N,as) = I^Cf^^ = Cg„{N,as) . (6) 

The matrix £ [N, Q/x, I^f) in Eq. (4) represents the evolution of the parton densities from scale 
io scale Q/x to NLL accuracy in In^ [23]. In effect, the combinations of parton distributions, 
evolution factors and matching coefficients follow the evolution of the initial-state hadrons from 
the QCD scale up to the scale at which the transforms suppress radiation. 

As indicated in Eq. (4), we adopt the renormalization scale Q/x in the C-cocfRcicnts. This 
is guided by standard Qr-resummation, where it was argued [12] that the scale should be pro- 
portional to 1/6. The choice of the scale in the C-coefficients actually only affects NNLL terms 
in the resummation. 

We note that in the inclusive limit 6 = the coefficients Cgjj and H do not quite reduce to 
the corresponding factor in threshold resummation. First, there is no analog of the nondiagonal 
Cg/q and Cgjq terms in the pure threshold case, since they are not associated with logarithmic 
behavior at large N. Second, the scale of ag in Cg/j reduces to Q/N rather than Q for 6 = 0. 
Finally, for the diagonal coefficient in our formulas above the combined effect at 0{as{Q'^)) is 
H X Cy/g X Cg/g ^ 1 H" ^ (n + SvT^), whcrcas the corresponding coefficient in the threshold case 
is [14] 1 + ^ (11 + 47r^). Such differences are present as well in the case of Z boson production 
through quark annihilation [23]. A development of the joint formalism to NNLL in both the 
threshold and small-Qr limits would account for them. 

We recall that organizing the C functions in the form of Eq. (4) in Mellin-A'" moment space 
has a great advantage: it enables us to explicitly separate the evolution of the parton densities 
between the scales fip and Q/x, embodied by the matrix S [23]. In this way, we can avoid the 
problem normally faced in Qt resummation that one needs to call the parton densities at scales 
far below their range of validity, so that some sort of "freezing" (or related prescription) for 
handling them is required. Our matrix £, on the other hand, may be expanded to LL and NLL 
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accuracy, consistent with all our approximations, as shown in ref. [23]. We then only need the 
parton distribution functions at the "large" scale yU^? ~ Q, whereas normally in resummation 
the product J2k^jk {N, l/h^iip) fk/H{N, /ip) is identified with fj/H{N, l/b). We also recall that 
the moment variable N and, as will be discussed below, also the impact parameter b are in 
general complex-valued in our approach, so that it is even more desirable to separate the complex 
scale Q/x from that in the parton densities. In practical applications, one usually has available 
only parton distributions in x-space. This obstacle can be overcome in various ways. We have 
described one practical approach in Appendix B of [23]. A very simple other method is to fit 
the parton distributions at scale /ii? in the x-region of relevance by a simple functional form that 
allows to take Mellin- moments analytically [28] . We have used both methods and found excellent 
agreement between the corresponding numerical results. 

The exponential E^^{N,b,Q, fi) in Eq. (1) is the Sudakov factor for Higgs production in 
joint resummation. It summarizes QCD dynamics in the transverse momentum range, Q/x ^ 
kr ^ Q, where real radiation is kinematically constrained and virtual corrections produce an 
exponentiating double-logarithmic suppression. We may think of it as interpolating between the 
scale of the partonic coefficients C and the short-distance annihilation process. E^J{N, b.Q, fi) is 
equivalent, with trivial changes in color factors, to the corresponding exponent in Z production 
[23], and is given by 



J5;/(iV.(-.«,,.) = -/J^^, 



(7) 



accurate to NLL level. Dependence on the renormalization scale /i is implicit in Eq. (7) through 
the expansion of aslkr) in powers of cks{ij) and logarithms of kT/lJ-- 

The anomalous dimensions in Eq. (7) are perturbative series in ag, Ag{as) — YliZo (^) -^g'\ 
and similarly for Bg . By examining the fixed-order cross section at NLO [2] the following familiar 
values necessary for NLL resummation can be derived [18, 29], 



4" 



Am 



Ca , 

Ca 
2 




4:TrNf) 



(8) 



where Ca = 3, Cp = 4/3, Tr = 1/2, and Np is the number of flavors. NLL accuracy requires 
using ^W, 5« and Af^ in Eq. (7). 

From an NNLO [3] calculation one can also determine a second-order coefficient for the 
function B{as)- This coefficient contributes to the cross section at the level of NNLL. Given 
our choices for H^^^ and hence for C^^\ and assuming the NLL expression for the Sudakov 
exponent (7), this coefficient is uniquely determined to be 



5f 



4 
3 



11^2 _ 3 
36 2^^ 



+ ^CpTrNp + CaNpTr {^-\ 



(9) 
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Even though the contributions associated with the coefficients H'^^\ C^^^ and B^^^ are all only 
of NNLL order, we will include them in our study. It was found in previous studies on Qt 
resummation for Higgs production that the combined contribution of these coefficients is actually 
numerically rather significant, partly due to the size of Ca relative to Cf- 

As was shown in ref. [19], there is actually some freedom in how one apportions contributions 
to H^^\ C^-^^ and B^^^ individually. A choice for one coefficient uniquely determines the others, 
but it is possible to shift terms among them, resulting in different but formally equivalent schemes 
for the resummation. We believe that our choice in Eq. (2) is well motivated physically, since 
with this choice the function H contains all genuinely hard short-distance effects at scales above 
Q. The function B^^\ on the other hand, enters with dependence on scales much smaller than 
Q, as can be seen from Eq. (7), so excluding hard contributions from B^^^ seems natural. We 
have nevertheless also investigated an alternative choice, for which all terms in H are absorbed 
into the C functions, so that H = 1. with B'^p adapted accordingly. This scheme has been more 
customary in previous studies. We found that our numerical results presented below do not 
change significantly if this "resummation scheme" [19] is used. We note that for our choice of 
H^^\ unlike for the H — 1 scheme, the corresponding numerical value for B^^^ is very small. A 
more detailed discussion of the interrelations between H^^\ C^^^ and B^^^ may be found in [19]. 

There are other NNLL contributions in the exponent, besides those associated with H^^\ 
C^^^ and B^^\ The third-order coefficient A^^^ also contributes at NNLL. Its value is so far 
known numerically [30], with its contributions proportional to Np and Np determined analyti- 
cally [31, 32, 33]. We have checked that when we include A^g^ at the level of its currently known 
approximation the impact on our numerical results is insignificant. 

We also note that in the context of joint resummation the form of Eq. (7) for the Sudakov 
exponent is itself only valid to NLL. The original jointly resummed Sudakov exponent given 
in [22] is, in fact, far more general than NNL and will correctly reproduce logarithms of any 
order for threshold and Qt resummation in the limits of large N and b, respectively. For large b 
and fixed N, the structure in (7) is correct to all logarithmic order, as is known from standard 
Qt resummation, where x ~ 6. In the opposite (threshold) limit of N large and b finite, we 
should in general introduce an additional function that reflects coherent soft radiation at wide 
angles, called in [10] and D in [33, 34]. For the purposes of this study, we follow standard Qt 
resummation as closely as possible and therefore use the form (7) also at NNLL. We emphasize 
at this point that a full investigation of NNLL effects in the joint resummation formalism would 
be very desirable. Our taking into account of H^^\ C^^^ and B''^ can therefore be only an 
approximation for the full NNLL effects. We defer a full NNLL analysis to a future publication. 

In the next subsection we discuss the Sudakov exponent in more detail. 
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2.3. The Sudakov exponent and the choice of x 



The exponent E^J{N, b, Q, ji) in Eq. (7) isolates the Sudakov suppression in the logarithm of the 
full NLL eikonal cross section in N- and 6-space [22], 



Jo (bkr) Ko I I + In ' 



Q 



Q 



-2 In AT J^^ '^A.iasikr)) , 



dkl 



(10) 



where here and below, we denote 



(11) 



with 7e the Euler constant. As in Ref. [23], we rewrite Eq. (10) in a form that is accurate to 
next-to-leading logarithm in both b and N, 



Ef^{N,b,Q,pi,piF) 



dkl 



+ 



QVx^ dkl 



Ag{as{kT)) hi\^j+Bg{as{kT)) 
- 2Ag(as{kT)) In (n) - Bg(as(kT))] . (12) 



Here we have added and subtracted a term with Bg, the term constant in N in the DGLAP 
splitting functions. To NLL, Bg coincides with the function defined in the previous section. We 
have also neglected logarithms of the factorization scale divided by Q. The Sudakov exponent is 
the first term in this expression. We will reinterpret the second term as the large- limit of the 
evolution of parton distributions in the next section. 

The function x(6(5, N) in Eq. (10) organizes logarithms of N and b in joint resummation [23] 
and is given by 

N 

x{bQ,N)^b + :^^ , (13) 



1 + ^^ 



where we define 



b = bQe'^/2 



(14) 



This form of x(iV, b) is chosen to reproduce threshold logarithms at NLL for large N at fixed 
b, and impact parameter logarithms at large b with N fixed. The functional form ensures that 
corrections at low Qt (hence, at large b) avoid unphysical singularities suppressed by 1/Qt 
compared to leading terms: as discussed in [23], any nonzero value for the parameter r] in Eq. 
(13) eliminates corrections of order in the large-6 limit, relative to leading behavior. It 
therefore does not introduce spurious logarithmic singularities in Qxda/dQ^dQl as Qt — > 0. To 
NLL, the cross section is independent of rj, as is the Qr-integrated cross section, found by setting 
b to zero. The shape of the cross section depends to some extent on the precise value of rj that 
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we choose, through subleading terms. For very large values of rj, x approaches its value in "pure 
recoil" 6-space resummation, x = b, although the integrated cross section is then ill-defined. 

In the following, we choose the value of ?7 = 1. We motivate this choice by comparing the 
expansions of the logarithm of x in terms of In b at fixed N, and of In N at fixed b. For rj — 1 
these expansions give rise to "power-hke" corrections of the form {N /Vf in the former case, and 
(h/NY in the latter, with no linear b/N corrections. We readily verify that the full NLL exponent 
(10) also has no corrections linear in h/N for h/N 0. The choice r] = 1 treats the remaining 
two sets of corrections, (N/b)^ and (6/7V)^, symmetrically. To estimate the sensitivity to rj, we 
will also exhibit results for = 1/2 and rj — 2 below. In Ref. [23], on the other hand, where we 
used the framework of joint resummation to analyse Z production at the Tevatron, a choice of 
T) = 1/4 was made. It led to a relatively simple structure of singularities in the complex b plane. 
We checked that the 1] dependence of the Qr distribution for Z bosons was insignificant [23]. 
This is, however, not necessarily the case for Higgs production at the LHC. The treatment of 
the inverse transforms (i.e. contours in the complex N and b spaces) for differing values of rj is 
discussed in section 3. 

The expansion of EjJ^ into leading, next-to-leading, and next-to-next-to-leading logarithms 
gives the form [22, 23] 



^PT,NNLL^^^ ^) ^ __ ^(0)^^) + 3 /i^ (/5, Q, /.) + 2 «,(/.) /if (/3, Q, /.) , (15) 

where each function h^g^ is a power series in 

P = boas{n)ln{x) , (16) 

with 

11Ca-4TrNf 
= 12^^ ■ ^'^^ 

The leading- and next-to-leading logarithmic functions /if and /if have already been presented 
in [22, 23]. For convenience, we provide them again in Appendix A. As stated above, we do 
not attempt a full NNLL study of the resummed cross section, but we do wish to include those 
NNLL terms that are known to produce numerically significant effects. At the level of the 
Sudakov exponent, this term is given by the contribution involving B^^^ and gives rise to the 
function /if: 

where the ellipses denote further, neglected, terms proportional to Af ,i?f ,Af , as well as to 
the third-order coefficient ^f . As mentioned above, we have checked that those terms among 
the neglected terms that are fully known lead to numerically insignificant effects. 
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2.4. Relation to threshold resummation 



In the resummed cross section, Eq. (1), the partonic functions Cg/HiQ,b, N, /i, /ip), Eq. (4), in- 
corporate the parton distribution functions fj/H{N, fip), the evolution matrices £jk {N, Q/x^ f^p), 
and the standard C coefficients coming from matching the fixed order to the resummed result. 

Within a strict threshold approximation, only parton-diagonal contributions to evolution are 
taken into account, which for Higgs production means the sole use of the large- A?" hmit of the gg 
anomalous dimension. The matrix £ {N, Q/XiI^f) then only has an entry in the gg position and 
reads _ _ 

^ ' 



^{N, QIx^[If) = exp 



-2AWlniV-5W , ' 



27r6o 



1 



where [23] 



= In 



a 



ln(l - 2/3) + NLL . 



(19) 



(20) 



In Eq. (19) we have for the sake of brevity only written down the LO/LL part of the evolution; 
extension to NLO/NLL is straightforward and included in our analysis. The £ in Eq. (19) may 
then be combined with the Sudakov exponent in (15) at /3 = ftoctsln -/V (i.e., 6 = 0) to reproduce 
the familiar exponent of Higgs threshold resummation [6, 13, 14]. 

However, as we discussed in [23], one can improve the above near-threshold approximation 
and use the full NLO singlet evolution matrix in £, with full A^-dependence for the anomalous 
dimensions. This amounts to the replacement 



-2Ag{as)\n(N) - Bg{as) — > 7iv(a.) , 
or equivalently, in terms of parton distribution functions. 



(21) 



5jg exp 



-2^W In AT 



27r6o 



5(1) 



fg,H{N, lip) (TV, Q/x, /.^) U,h{N, ixf) . (22) 



In this way, £, and hence the functions C in Eq. (1), obtains a full matrix structure and becomes 
an ordered exponential. On the other hand, the standard threshold resummation formula is 
recovered from Eq. (1) by integrating over Qt to set 6 = 0, and then undoing the replacements 
of Eqs. (21) and (22). 

One of the assets of the joint resummation of Eq. (1), incorporating (22), is to resum coUinear, 
but non-soft, terms of the form a" ln^"~^(x)/iV to all orders [23]. Such terms are obviously 
subleading with respect to the threshold approximation, but may have some relevance for phe- 
nomenology. They were discussed in [6, 13, 14] in the context of threshold resummation, where 
they occur as a^\ri'^~^{N) / N . In this formalism they are brought under control automatically. 
In our analysis of Higgs production, we actually find that other formally subleading terms can 
also be rather important. We will discuss this issue in section 5 below. 
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3. Inverse transforms and matching 



We have already mentioned in the introduction the importance of suitably defining the inverse 
Mellin and Fourier transforms. We follow here closely our previous phenomenological study for 
elcctrowcak annihilation [23], where all details may be found. Let us briefiy recall the main 
ingredients. 

The contour for the inverse Mellin transform is chosen according to the "minimal prescription" 
of [25]. It is parameterized as 

AT = C + ze^'^ , (23) 

where the upper (lower) sign applies to the upper (lower) branch of the contour, with < z < oo 
(oo > 2; > 0). The constant C is required to lie to the right of the rightmost singularity of the 
parton distribution functions. At any finite order in perturbation theory, all values of C > 0, and 
all angles tt > > 7r/2 are equivalent. In the resummed cross section, however, the singularity 
s-t X = Pl = exp[l/26otts(Ai)] introduces an ambiguity in the transform, which may be resolved 
by choosing C < pl [25]. 

The inverse Fourier integral is also defined as a contour integral [23, 26]. It is performed using 
the identity 

/"OO roo 

d%e"^-'' f{b) = 2tt dbb Jo{bq) f{b) = 71 dbb [hi{bq,v) + h2{bq,v)] f{b) , (24) 
Jo Jo 

and employing Cauchy's theorem to deform the integration over real b into a contour in the 
complex b plane [26, 23], as shown in Fig. 1. Here the auxihary functions hi^2 [26] are related 
to Hankel functions. They distinguish between the positive and negative phases in Eq. (24). 
The b integral can thus be written as a sum of two contour integrals, of the integrand with hi 
(/i2) along a contour in the upper (lower) half of the b plane. The precise form of the contours 
becomes unimportant as long as the contours do not run into the Landau pole or singularities 
associated with the particular form (13) of the function x- One also needs to make sure that 
the lines of singularities do not cross, which would make it impossible to draw a fixed contour 
without running into a line. This can be ensured by choosing the opening angle of the N contour 
to fulfill the condition 7i/2 < (p < 71 - tan"^ {y/Hv - 1|). 

This choice of contours in complex transform space is completely equivalent to the original 
form, Eq. (24) when the exponent is evaluated to finite order in perturbation theory. In the pres- 
ence of the Landau pole arising in the resummed formula, it is a natural extension of the iV-space 
contour redefinition above [25], using a generalized "minimal" exponent (15). We emphasize that 
joint resummation with its contour integration method provides an alternative to the standard 
b space resummation. Joint resummation has built-in a pcrturbativc treatment of large b values, 
eliminating the need for a 6^, [12] prescription for the exponent, or for a freezing of the scale 
of parton distributions at large b or low Qt- In this way, we can derive entirely perturbative 
cross sections for Higgs production. We note that this definition of the 6-space contour was also 
adopted in the recent study [17], while in [27] an extrapolation of the perturbative exponent to 
large b is used to define the inverse transform. 
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Figure 1 : Choice of contour for b integration (thick sohd hues) for rj = 1. The straight sections of 
the contour from to be are to be interpreted as on the positive real axis. The remaining curves 
represent hues of singularity discussed in [23] . 

At large values Qt ~ rrih, wc need to match the resummed cross section to fixed-order 
perturbation theory. This is achieved in the following simple way [17, 23]: 

da _ da'"' rfa^^P^'^) da^"""^^^^ 
dQMQ^T ^ dQ^dQl ~ dQ'^dQl ^ dQ^dQl ' 

where da'^/dQ'^dQ'^ is given in Eq. (1) and da^^^^^^dQ'^dQ'^ denotes the terms resulting from 
the expansion of the resummed expression in powers of Q;s(/i) up to the order k at which the 
fixed-order cross section da^^^'^^^^ dQ'^dQ'^ is taken. The last two terms on the right-hand-side 
of Eq. (25) should cancel each other at small Qt, since order by order the terms generated by 
the resummed formula must reproduce the singular behavior of the cross section at small Qt- In 
our study, we do the matching using k = 1 (see [7]) in Eq. (25). 

4. Predictions for the resummed Higgs cross section 

For our calculations of the Higgs transverse momentum distribution at the LHC, obtained in 
the framework of joint resummation, we choose rrih = 125 GcV and the factorization and renor- 
mahzation scales ji — jip — Q = rrih. We use the CTEQ5M [35] parton distribution function 
parameterizations. The contour parameters are = 17/327r and (pb — (f) + tan~^ ^\4:ri — 1\, with 

77 = 1. 

Fig. 2 shows our result for da/dQT as a function of Qt- The solid curve shows the result 
obtained for the formulas given in Sec. 2. For comparison, we also display the result for "pure- 
Qt" resummation, i.e., for the case when we do not take into account threshold resummation 
effects, X = b (dashed line). We find that the corresponding result is somewhat higher and 
broader than our jointly resummed curve. We can compare this curve to previous analyses in 
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the literature [9, 15, 16]. It turns out that our "pure-Qy" resummed result is about 30% higher 
(in the peak of the distribution) than the one of the Resbos [15] code, and about 20% higher 
than the result of [16]. On the other hand, our distribution is lower by about 15% than the NLL 
result of [17], matched with the NLO predictions. These differences appear to be due in part to 
different treatments of the large-6 region in the various formalisms, and to different organizations 
of the coefficients. 




Figure 2: Transverse momentum distribution for Higgs production at the LHC in the framework 
of joint resummation and of "pure-Qr" resummation. 



The jointly resummed distribution in Fig. 2 was obtained after adding a small non-perturbative 
term of the form —g h"^ to the exponent. As discussed in [23], it is an advantage of joint resumma- 
tion that no non-perturbative contribution is required in order to obtain results for any nonzero 
Qt- The reason for introducing a non-perturbative component here is purely technical. Our 
choice of ?7 = 1 necessitates small values of angle 0, making it time consuming to achieve numeri- 
cal stability at low Qt- The value of the parameter g was adopted from the study [36], where an 
analysis of Upsilon hadropro duct ion was performed as a means of studying the non-perturbative 
contributions in processes with two gluons in the initial state. We use the value g = 1.67 GeV^ 
corresponding to the lower end of a range of values extracted from data. In any case, on the 
basis of previous studies of the dependence of 6-space resummation for Higgs production on the 
non-perturbative contribution [16, 36, 37], we expect the dependence on the parameter g to be 
weak for values of Qt around or above the peak of the cross section. In effect, a nonzero value 
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of g simply enables us to extrapolate the cross section to Qt = 0. 

In section 2 we introduced the parameter rj in the definition of the function x, Eq. (13), and 
as discussed above chose 77 = 1. To investigate the dependence of the results on the value of this 
parameter, we computed the Qt distributions for 77 = 1/2 and rj = 2. As can be seen from Fig. 3 
the numerical sensitivity of the results to the value of rj is rather small in this range. 




20 40 60 80 100 120 

Qt (GeV) 



Figure 3: Dependence of the jointly resummed transverse momentum distribution for Higgs 
production at the LHC on the value of the parameter 77. 

5. Joint resummation and the role of non-threshold 
effects 

In a sense. Fig. 2 completes our study, since it predicts a cross section that remains to be 
measured. Nevertheless, a number of cross checks in joint resummation suggest that Fig. 2 is not 
the final word, and that to get a full picture of the Higgs cross section in Qt it may be desirable 
to extend the formalism to include the contributions that give rise to dominant small-x behavior. 

In this section, we compare the jointly resummed cross section in more detail to existing 
fixed-order and threshold-resummed calculations. 
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5.1. NLO expansion in Qt space 

Fig. 4 shows the difference da^^^'^^^'^ / (IQt — do^^"^^^^ / cLQt and the fractional deviation 



A 



^^fixed{l) rf^exp{l) 



dQ 



T 



dQj 



I 



^^fixed{l) 

dQr 



(26) 



using nih = 125 GeV, scales ^ = fip = Q = fnt, and the CTEQ5M [35] parton distribution 
functions. We find very good agreement in the region Qt < 10 GeV where resummation is 
necessary. Beyond 10 GeV, the agreement is naturally less exact but still reasonable. This is the 
region where matching to finite order is appropriate, as expressed by Eq. (25). From the figure 
we see that the 'recoil' approximation, x — > 6, is slightly more accurate at low Qt, while the 
joint form remains stable at larger values of transverse momentum. The somewhat lower value of 
the expanded joint cross section is due in part to the lower scale at which the gluon distribution 
is evaluated compared to pure Q^-resummation {Q/x as opposed to l/b). At larger values, the 
stability of the jointly resummed cross section makes it easier to match to fixed order (a property 
shared with the Q^-resummed treatment in Ref. [17]). 




Qj (GeV) 



Figure 4: Fractional deviation A (as defined in Eq. (26)) between the "exact" 0{as) result and 
the 0{as) expansion of the jointly and the pure-Qr resummed cross sections. 



5.2. Total cross section and subleading terms in N space 

When integrated over Qt, the jointly resummed distribution formally reduces to the threshold 
resummed cross section in the large- limit. When we numerically compare the two formalisms. 
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however, we find that the total cross section in joint resummation is around ten percent lower 
than pure threshold resummation at NLL: after matching to NLO we find the NLL threshold- 
resummed cross section to be 39.4 pb (consistent with results shown in [14]), whereas our jointly 
resummed cross section integrates to 35.0 pb. This slight disagreement implies that terms that 
are subleading with regard to the threshold approximation play a fairly important role. Since 
our formalism has been designed to automatically reproduce all behavior singular at Qt 0, 
we are led to associate the "missing area" with contributions to the Higgs cross section at rather 
high Qt. 

It is instructive to investigate these points in more detail. To this end, we expand the jointly 
resummed, and Qr-integrated, cross section to fixed order as (not counting the overall of the 

cross section), and compare it to the "exact" expressions for the Higgs production cross section at 
this order. For simplicity, wc set the renormalization and factorization scales to m/j. Expansion 
of Eq. (1) to 0{as) gives for the partonic cross sections in the gg and gq scattering channels: 



^99 



(jf — |-4Ca In^ N + 87r6o In AT + 11 + 37r^ 

ACa , 4C^ 



-2 In AT 



+ 



r(0) 



etc 



N{N-1) (iV + l)(iV + 2) 
+ AT + 2 1 



4CA^i(Ar)+47r6o 



''^"'2^^^r^^''%(7V + l)(Ar-l) ■ AT + l 



+ 



(27) 



where 5i(iV) = Ej=iJ~^ = V^(A^ + 1) +7b, with ^ the di gamma function. Since we are interested 
in the near-threshold region, we can expand this further to the large- iV limit. Using ip[N -|- 1) = 
lnA^+ l/(2iV) + 0(l/iV2), we find: 



\99 



^91 



a.. 



27r 



4CAln^A^ + 4CA 



(0) 



27r 



N 




+ 



(28) 



The full ( "exact" ) NLO cross total section has been calculated in refs. [2] . When taking MeUin-A/" 
moments of these results, we obtain 



-99 
^ exact 



exact 



a^^^ ^ UCaSUN) - 8Ca 



1 



1 



3C, 



+ 



13Ca 



A^(A^ + 1) (A^-l)(A^ + 2) 



N{N-1) {N + l){N + 2) 



Si{N) + 11 + 47r^ 



1 1 

+ 



^ 27r 



+ 



TV^ TV + 2 



(AT -1)2 7V2 (AT + 1)2 (A^ + 2)2 
2 2 1 



N{N + 1){N -1) (AT -1)2 (AT + 1)2 

1 1 1 



2(Ar-l) N 2(Ar+l) 



(29) 
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At large N, this gives 



»" = <'Sc4[-21„^+l]+o(l^) . (30) 

Comparing Eqs. (28) and (30), we see that they agree even to the subleading order 0{1/N), 
except for the contribution proportional to tt^ discussed in section 2.2. As noted there, had we 
fully developed the joint formahsm to NNLL for both the threshold and the small-Qr limits, 
it would automatically account for this difference in the C coefficients. At the level of our 
rudimentary inclusion of NNLL effects at large b, we need to accept the mismatch in the tt^ term. 

In addition, there are differences between the 0{as) expansion of the resummcd cross section, 
Eq. (27), and the exact 0{as), Eq. (29), that are formally even more suppressed than the 
above tt^ term, or even the 0{laN/N) terms mentioned earlier, but that still have a fairly 
large impact on the results. These correspond to the terms oc 1/(A^ — 1)*^ {k = 1, 2) in 
these equations. For instance, there is a term —8Ca^^N/N{N — 1) in Eq. (27) that arises 
from our treatment of evolution. In Eq. (29), on the other hand, we find the combination 
-SCaSi{N)/N{N - 1) + 13Ca/{N -1){N + 2) + ACa/{N - if. Numerically, the former term 
leads to a sizeable negative contribution in Eq. (27), while the latter terms cancel each other 
to a large extent in Eq. (29). This is the main reason why our jointly resummed cross section 
integrates to a smaller area than one obtains for the standard threshold-resummed cross section 
on the basis of Eq. (30), where of course no oc 1/{N — l)'^ contributions are present since they 
are suppressed from the threshold point of view. 

The relevance of the 1/(A'" — l)'^ contributions results from the fact that these terms actually 
constitute the most singular behavior in the limit opposite to partonic threshold, z = Q'^/s 0. 
For Higgs masses of a few hundred GeV and LHC's center-of-mass energy one is actually quite 
far from threshold. We note that at higher orders in ag the leading terms singular at = 1 
generalize to ag/(N — 1)^*^, corresponding to ag\n'^'^~^{z)/z in z space. These terms, too, can 
be resummed to all orders in ag, and the resummation was achieved in Ref. [38]. To what 
extent it would be possible to incorporate this small- 2; resummation into our joint formalism is 
an interesting question that, however, is beyond the scope of this paper. That the leading and 
subleading terms oc 1/{N — 1)^ {k = 1, 2) almost entirely cancel numerically in the full NLO 
inclusive cross section, Eq. (29), is rather remarkable, albeit probably not accidental. It is only 
because of this feature that the near-threshold approximation for the NLO cross section, Eq. (30), 
becomes a good approximation. 

Assuming that threshold resummation indeed does yield a reliable answer for the integrated 
total cross section, it becomes then an interesting question what the implications for the Qt- 
dependence of the Higgs cross section will be, if it is derived in such a way that it the area 
underneath the distribution integrates to the threshold-resummed cross section. A possible, but 
clearly not unique, way of making up for the "missing" area is to redefine the C coefficients of 
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Eq. (5) in the following way: 




13C, 



Cgig {N, as 



+ 



AC, 



(N-l){N + 2) (iV-1) 



(31) 



Adding terms of this form to the C coefficients is completely consistent with resummation for 
both the threshold and Qt resummation cases. The term 1/(1 + 6^) is designed so that it becomes 
unity in the Qr-integrated cross section and does not produce behavior singular in Qt- We do 
not claim that the coefficients in (31) are necessarily more realistic or better motivated than the 
standard ones of Eq. (5); we rather regard them as one model for the effects associated with 
the cx 1/{N — 1)^ terms on the Higgs Qt distribution, for which the QT-integrated cross section 
becomes numerically equivalent to a NLL threshold-resummed cross section. 

In Fig. 5 we show the effect of the modified C coefficients on the Higgs Qr-distribution. 
As expected, there is no difference between the Qt distributions with the C coefficients in the 
forms (5) and (31) in the region where the distribution peaks; however, the new coefficients cause 
the Qt distribution to be higher at intermediate / large values of Qt- The value of the total 
cross section is now about 41 pb. 



6. Conclusions 

In this paper we have presented a study of the transverse momentum distribution of Higgs bosons 
in the context of the joint resummation formalism. Guided by the literature on the more standard 
Qt resummation, we have further developed the joint resummation formalism to make it suitable 
for application to the Higgs production cross section. This included the implementation of some 
NNLL terms that have a significant numerical effect on the Higgs Qt distribution. 

We have found in our study that effects associated with threshold resummation have a modest 
importance for the Higgs cross section at small to moderate Qt- For this region, where the 
distribution peaks, our results further confirm the applicability of standard Qt resummation. 

On the other hand, threshold resummation is very relevant for the Qr-integrated cross section, 
as is known from previous studies in the literature. For joint resummation, the integrated cross 
section is formally identical to the threshold resummed one in the large- limit, but contains 
non-leading contributions that go beyond standard threshold resummation and which, as we 
found, decrease the cross section. They are associated with small-x terms from parton evolution, 
in a sense of opposite origin to the threshold logarithms. We have argued that their presence 
indicates that the Higgs Qt distribution at rather large Qt receives sizeable contributions also 
from sources not related to Sudakov logarithms. Our study also suggests that it will be useful to 
extend joint resummation to NNLL. Such an extension will help clarify technical issues, such as 
the choice of scales in coefficients and parton distributions, which have a modest but significant 
effect on the magnitude of the cross section even near the peak. 
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Figure 5: Effect of introducing extra terms in tfie C coefficients, as shown in Eq. (31) and 
discussed in the text, on the transverse momentum distribution for Higgs production at the LHC 
in the framework of joint resummation. 
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A Leading- and next-to-leading logarithmic expansion 

Here we give the functions h^g^ and h^^^ that determine the jointly resummed exponent (15) to 
leading and next-to-leading logarithmic accuracy [22, 23]: 

hfHP) = ^ [2/3 + ln(l - 2/5)] , (32) 
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2 ^ 1-2/5 



+ 



BO, 
27r6o 



ln(l - 2/3) 



+ 



27r6o 



^(2) 

7r6o 



2/3 



1-2/3 



+ ln(l - 2/3) 



(33) 



with the coefficients Bf from Eq. (8), /3 from Eq. (16), and 6o a-s in Eq. (17). In addition, 



17C1 - lOC^TRiV^ - QCfTrNf 
247r2 



(34) 
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